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ABSTRACT 

With the use of a detailed Milky Way nonaxisymmetric potential, observationally 
and dynamically constrained, the effects of the bar and the spiral arms in the Galaxy 
are studied in the disc and in the stellar halo. Especially the trapping of stars in the disc 
and Galactic halo by resonances on the Galactic plane, induced by the Galactic bar, 
has been analysed in detail. To this purpose, a new method is presented to delineate 
the trapping regions using empirical diagrams of some orbital properties obtained 
in the Galactic potential. In these diagrams we plot in the inertial Galactic frame 
a characteristic orbital energy versus a characteristic orbital angular momentum, or 
versus the orbital Jacobi constant in the reference frame of the bar, when this is the 
only nonaxisymmetric component in the Galactic potential. With these diagrams some 
trapping regions are obtained in the disc and halo using a sample of disc stars and 
halo stars in the solar neighbourhood. We compute several families of periodic orbits 
on the Galactic plane, some associated with this resonant trapping. In particular, we 
find that the trapping effect of these resonances on the Galactic plane can extend 
several kpc from this plane, trapping stars in the Galactic halo. The purpose of our 
analysis is to investigate if the trapping regions contain some known moving groups 
in our Galaxy. We have applied our method to the Kapteyn group, a moving group 
in the halo, and we have found that this group appears not to be associated with a 
particular resonance on the Galactic plane. 

Key words: Galaxy: disc - Galaxy: kinematics and dynamics - Galaxy: structure - 
galaxies: spiral structure, bars 


1 INTRODUCTION 
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them have emerged, from simple disruptions of open 
clusters by Galactic tidal effects, to cosmological ori- 
Jns, such as dwarf sat ellite galaxies or globular clusters 
[Penarrubia et al.l l2009l f . to perturbations by re sonances 
with the Galacti c bar and/or sp i ral a r ms ([Pehnen 
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together contribute to form the multiple moving groups 
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seen today in the Galaxy, but is there one more effective 
than the others, should we be seeing more groups produced 
by intrinsic secular processes than the ones produced 
by disrupted satellites or clusters, or can we identify 
clearly what phenomenon produces which particular group? 
Hopefully soon, these and many other questions will be 
better answered with the impending arrival of the great 
new large surveys (such as Gaia, SLOAN/APOGEE, the 
Geneva-Copenhagen Survey, RAVE, SEGUE, etc.) joined 
to the best theoretical tools to elucidate the different 
processes. 


Meanwhile, one of the intricate problems to solve today 
is the detection of these groups. Some of the best known 
methods to identify moving groups in the literature are those 
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2 Autores 


S Also the method of Lvnden-Bell Lvnden-B^ 
and of iLvnden-Bell I (|l999l l that analyses positions 
and velocities of globular clusters and satellite galaxies, 
searching for halo fossil tidal streams that mark orbits 
followed by satellites merging with the Galaxy. Another 
me thod is the one known as Great Gircle Gell Gounts (GG3) 
by I Johnston et al.l (|l996h , that makes use of the fact that 
orbits dominated by the spherical outer halo conserve ap¬ 
proximately their orbital plane orientation and leave behind 
‘great circles’ on the sky. And, methods to search for mov¬ 
ing groups in the disc, that work in velocity space, con¬ 
sidering that all the stars in a given moving group move 
toward the same direction fe.g. iHoogerwerf AguilanI 19991 : 

Ide Bruiinn1ll999l b A method based on t he lumpiness in in¬ 
tegral s of motion has been employed bv iHelmi de Zeeuwl 
tood ) to study the merger history of our Galaxy. For a re¬ 
view of methods, see this last reference. 

Previous work has demonstrated the influence of galac¬ 
tic bars on resonant trapping of orbits on the Galactic 
plane and out of it. For example particles on bar reso¬ 
nances can exchange or transfer angular momentum to par¬ 
ticles in the dark matter ha l o, dr i ving dynamical evolu¬ 
tion f (iLvnden-Bell KalnaisI 19721: iTremaine Weinberg 


19841: lAthanassoulal l2002l: iMartin ez-yal Duesta et al.l l200d 


Sellwood I |2014 iGeverino Klvpinl 2QQ7l l and references 
therein). Although N-body simulations will probably be 
the best test-bed models to produce studies like the one 
presented here, it is still intricate to extract information 
f rom resonant particles b ecause of their lack of resolution 
f (|Weinberg Kat3l2007c] li3B. Therefore, we employ in this 
work a very detailed Milky-Way steady model that is de¬ 
scribed in Section [3l 

Assuming that a large fraction of moving groups can be 
explained by intrinsic properties of the Milky Way Galaxy, 
such as its nonaxisymmetric potential structure, in this work 
we analyse the trapping of stars by resonances on the Galac¬ 
tic plane produced by these nonaxisymmetric components, 
and our aim in this and in a second paper is to investigate the 
possible relation of these trapping regions with some known 
moving groups in our Galaxy. For this purpose, we have de¬ 
veloped a new straightforward tool to analyse the resonant 
trapping. The method requires the computation of stellar 
orbits in a nonaxisymmetric Galactic potential, and makes 
use of diagrams of a ‘characteristic’ orbital energy versus a 
‘characteristic’ orbital angular momentum, or versus the or¬ 
bital Jacobi constant if the Galactic potential employs only 
the bar in its nonaxisymmetric components. 

A sample of halo stars and disc stars in the solar neigh¬ 
bourhood is employed to show how the proposed method 
works. Our Galactic potential uses mainly the bar as the 
nonaxisymmetric component, but results are given showing 
the additional effect of the spiral arms. An axisymmetric po¬ 
tential is also considered, to compare with the main features 
that arise in the nonaxisymmetric case. 

This paper is organised as follows. In Section [2l the 
observational stellar sample is described. In Section [3l the 
properties of the three-dimensional Galactic potential em¬ 
ployed to compute orbits are summarised. In Section |4] the 
new method and its technique to analyse resonant trapping 
is introduced. Section [5] gives a study of periodic orbits on 
the Galactic plane obtained in a Galactic potential including 
only the bar in its nonaxisymmetric components; this study 


gives the orbital support on which the proposed method is 
based. Section [6]shows the important influence that resonant 
families on the Galactic plane can have on stellar orbits at 
different ^-distances from this plane, which is fundamental 
for the application of the method to stars in the Galactic 
halo. In Section [7| this method is complemented by present¬ 
ing the spectral analysis of the Galactic orbits of the sam¬ 
ple stars and of the computed periodic orbits. In Section [8] 
the method is applied to the Kapteyn moving group in the 
Galactic halo. Finally, Section [9] presents our conclusions. 


2 THE OBSERVATIONAL SAMPLE OF STARS 

The sample used in this work consists of 1642 halo and disc 
stars from the catalogues of uvby-^ photometry of high- 
velocity and metal-poor stars observed at the National As¬ 
tronomical Observatory, San Pedro Martir, Mexico (here¬ 
after SPMl (IS chuster NissenI Il988l : ISchuster et al.l 1 19931 . 
l200fit ISilva et, al.ll2012h with 711. 553, 442, and 143 stars 
in each of these catalogues, respectively. However, not all 
of these stars have complete kinematic data, and some of 
these stars have been repeated in more than one catalogue, 
w hich has been ame nded in the hnal data base as explained 
in ISilva et ah (l2012h. Also, about ha lf of the stars in the first 
catalogue (jSchuster NissenI Il98^ i have been observed at 
the Danish telescopes. La Silla, Ghile, (hereafter LS): 36.0% 
of the stars in this first catalogue were observed exclusively 
in the north at SPM, those with the note ‘N’ in their Table 
VIII; 50.1% exclusively in the south at LS, those with ‘S’; 
and 13.8% at both observatories, denoted with ‘N/S’. 

The coordinates and proper motions employed i n this 
pape r have been taken primarily from: Hipparcos (|ESA I 
Il997li. Tycho-2 ( Hpg et all l200ol i , and the revised NLTT 
(|Salim &: Gould] 2003l i. Fot a few stars other sources have 
been u sed for the pr oper moti ons, such as th e original 
NLTT (|Luvten Ill979al i. the LHS (iLuvten lll979lJ), the Low - 
ell Proper Motion Survey (iGiclas e t al.lll9 59l. Il96ll. j^75l i. 
and the PPM Star Gatalogue ([Roser Bastian Il99lh ~ln 
general the proper-motion errors range from 0.5-3.0 mas/yr 
per component for the best sources, such as Hipparcos and 
Tycho-2, to as large as > 10 mas/yr per component for the 
other sources. The median errors of the proper motions are 
±1.51 mas/yr for /Xq;, and ±1.20 mas/yr for ys. 

The distances were derived f rom the phot o metric 
absol ute-m agnitude c a librat ion of ISchuster et al.l (|2004l . 
l200fih and ISilva et, all (12012^ . or from parallaxes in SIM- 
BAD for stars with errors less than 10% or for stars outside 
the colour range of the photometric calibration. The median 
error of these distances is ±6.8 pc, or ±9.6% of the distance. 

The radial velocities have been sel ected from the litera ¬ 
ture, from a number of so urces , su ch as Garnev et al. (Il994ll: 


iBarbier-Brossat et al. (Il994 b Barbier-Brossat Figonl 
(l2000h: Ryan Norris I (|l99lh ; iFouts Sandage I (|l986h : 


lAbt Biggs I (|l972l l and Nordstrom (private communi¬ 
cation). SIMBAD has been consulted to obtain the most 
precise and up-to-date values possible; about 20% of our 
radial velocities ha ve been updated from our previous study 
(ISilva et alJ 12012lV The errors of these radial velocities 
range from a few tenths of a km for the best sources 
(Garney et al. 1994; Nordstrom, private communication). 
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to ~ 7 km s ^ for the older sources. The median error of 
the radial velocities is ±1.45 km s~^. 


3 THE GALACTIC MODEL 

In order to construct a comprehensive orbital study of stars 
in ‘moving groups’ in the Galactic disc and halo, a detailed, 
observationally restricted, semi-analytic (i.e. with some of 
the functions numerically solved), three-dimensional model 
of the Milky Way galaxy, composed of axi symmetric and 
nona xisymmetric potentials, is employed (| Pichardo et al.l 
l2012l b In this model, some modif ications are made on th e 
axisymmetric Galactic potential of lAllen SantillanI (|l99lh . 
This axisymmetric potential consists of three components: 
a disc, a spherical bulge, and a massive spherical halo. The 
first modifications to build the nonaxisymmetric Galactic 
potential consist in that all the mass in the spherical bulge 
is employed to build the Galactic bar, and a small frac¬ 
tion of the total mass of the disc is employed to build the 
three-dimensional spiral arms. Therefore, only two axisym¬ 
metric components are conserved in the final model, the 
slightly diminished disc and the spherical halo. A boxy bar 
is employ ed to represent t he Ga lactic bar; its model is de¬ 
scribed in IPichardo et all (|2004l i. The model for the three- 
di mensional spiral arms, called PERLAS, is described in detail 
in I Pichardo et al.l (|2003l i . The Galactic potential has been 
rescaled to the Sun’s galactocentric distance, Ro = 8.3 kpc, 
and the local rot a tion v elocity, ©o = 239 km s“^, given by 
iBrunthaler et al.l (|20Ilh . 

A brief description of the observational/theoretical pa¬ 
rameters utilised to restrict the nonaxisymmetric compo¬ 
nents of the model is presented next and summarised in 
Table [1] (further details on the parameters of the mod el and 
restrictions were introduced in IPichardo et all (|20I2l P. 


3.1 The Galactic Bar 


In I Pichardo et al.l (|2004l l three different models for the 
Galactic bar were introduce d , whic h approximate the den¬ 
sity model of IPreudenreichl (|l998l ) of GOBE/DIRBE ob¬ 
servations of the Galactic centre: a prolate, a triaxial 
and a triaxial boxy bar. Eor the computations in this 
work, we have adopted the triaxial boxy bar model, that 
fits better the Galactic observations. Table [1] summarises 
the bar’s structural parameters, such as: the major semi¬ 
axis, scale lengths, and axial ratios. The considered to¬ 
tal mass for the bar, 1.6 x I0^° M©, is within obser¬ 
vational e stimations t hat lie in t he range I ~ 2 x I0^° 
M© ( e-g-. iKent 199d: Ztm3 1994 Dwek et~al]|l995 : iBlumI 


1 1 9951 : l^anek et al.l 1 19971 : I Weiner SellwoodI Il99f ). The 
angular speed of the bar and the present orientation of 
its major axis are still controversial fe.g.. [ Bi nnev at al.l 
199ll: IW einer Sellwoo dI Il999l: i Bissantz Ger hardI l2002l: 


Babusiaux Gilmord l200^ Minchev et ahl ' 20071 : iGerhardI 

2002l i . We consider 0 = 20° for the present-day orientation of 


the major axis of the Galactic bar; this is the angle between 
this major axis and the Sun-Galactic centre line. A long list 
of studies have been presented in the literature that estimate 
the most important dynamic al parameter of the Galactic 
bar, its angular velocity, (|Gerhardll201lL and references 
therein). iGerhardI (|201l| j concludes from his review that the 


most likely range in Qb is 50 - 60 kms“^kpc“^. Eor our 
computations, the value Qs = 55 kms“^kpc“^ is employed, 
as listed in Table [U Thi s value is consistent with the recent 
estimate of given bv lAntoia et al.l (|2014l b Qb = (1.89 ± 
0.08)flo, with flo the local angular velocity; in our potential 
this gives = 54.4 ± 2.3 kms“^kpc“^. 


3.2 The Spiral Arms 

The spiral-arm model employed in this work is the one 
called PERLAS (| Pichardo et al.l l2003l j: it is formed by a 
three-dimensional bisymmetric steady potential consist¬ 
ing of a superposition of inhomogeneous oblate spheroids 
along a logarithmic spiral locus, adjustable to better rep¬ 
resent the available observations of the Galactic spiral 
arms. This simulates the main Galactic spiral arm s as in 
iBeniamin et al.l (|2005l j and iGhurchwell et al] (|2009l j , based 
on the Spitzer/GLIMPSE database. The position of the 
spiral arms at the p resent time is taken as in figure I in 
I Pichardo et al.l (|20I2l b i.e. the line where these arms emerge 
at the inner Galactic region lags behind the major axis of 
the bar making an angle of 40° with this axis. The spe¬ 
cific parameters of the spiral arms, such as angular velocity, 
pitch angle, etc., are provided in Table [J an d more specihc 
detail s of the parameters are presented in I Pichardo et al.l 
(|20I2l i. The spiral-arm mass is distributed as an exponen¬ 
tial decline along the arms, and their strength is related 
with their total mass, which is a small fraction of the disc 
mass. To quantify the st r ength of the arms, the fun ction Qt 
(| Sanders L Tubbj II98QI : iGombes Sanded Il98lh is com¬ 
puted. The maximum value of Qt over the radial extent of 
the spiral arms, called Qs, is a me asure of the stren gth of 
the spiral arms. In agreement with iButa et al.l (|2005l j. and 
considering that for Sbc galaxies, as for our Galaxy, Qs is ap¬ 
proximately less than 0.25, the corresponding range in total 
mass for the spiral arms, calculated from PERLAS model, is 
obtained with Marms/Mdisc = 0.04±0.0I. We have taken the 
central value 3.9 x 10^ M© for the mass of the spiral arms. 
Einally, in the same manner as for the Galactic bar, there 
are different methods to determine the spiral-arm angular 
velocity; for a review see IGerhardI (|20Ilh . We use the mean 
value of Qs = 25 kms“^kpc“^, provided by these methods, 
as listed in Table [U 


4 A NEW METHOD TO ANALYSE 
RESONANT TRAPPING 

The orbits of the 1642 stars in the catalogues of Section [2] 
are computed for different forms of the Galactic potential, 
such as (I) an axisymm etric case (the rescaled potential of 
lAllen Santillai] (|l99lh b (2) a nonaxisymmetric case that 
includes the disc and halo as the axisymmetric components, 
plus the boxy bar, and (3) the nonaxisymmetric case (2) 
with the addition of the spiral arms to its nonaxisymmet¬ 
ric components. All the orbits are computed backward in 
time up to I0^° yr. Some parameters are obtained in each 
orbit, such as the orbital energy per unit mass, E, and the 
^-component of angular momentum per unit mass, h (the z 
axis being perpendicular to the Galactic plane); these two 
important parameters in our method are obtained in all the 
employed forms of the Galactic potential, and are computed 
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Table 1. Parameters of the Nonaxisymmetric Galactic Components, 
Sun’s Galactocentric Distance and Local Circular Rotation Speed 


Parameter 

Value 

References 

Triaxial Boxy Bar 

Major Semi-Axis 

3.13 kpc 

1 

Scale Lengths 

1.7, 0.64, 0.44 kpc 

1 

Axial Ratios 

0.64/1.7, t)AAll.7 


Mass 

1.6 X lO^o M© 


Angle Between Major 

Axis 


and the Sun-GG Line 20° 

2 

Pattern Speed (H^) 

55 kms ^kpe ^ 

3,7 

Spiral Arms 

Pitch Angle (z) 

15.5° 

4 

Scale Length (H^) 

3.9 kpc 

5 

-Marms/ iodise 

0.04 


Mass 

3.9 X 10® M0 


Pattern Speed (H^-) 

25 kms“'^kpc“'^ 

3 

Rq 

8.3 kpc 

6 

00 

239 km s“^ 

6 


(1) l^eudenreicyi 1998lh (2) iGerha xdl (l2002h: (3) iGerhardl l^2011^: 
l'4~l iDrimmell ll20nnh~ iBcnjamin et a] 1 (120051) : (6) iBrnnt haler et al l 
ll201iri: r74Antoia et alJ ll2014h. 


in the Galactic inertial frame. Other computed parameters 
are the Jacobi constant per unit mass, J, (in the case of the 
nonaxisymmetric potential including only the boxy bar in 
its nonaxisymmetric components; case 2 above), the mini¬ 
mum and maximum distances from the Galactic centre, the 
maximum distance from the Galactic plane, and the orbital 
eccentricity. 

Figure [T] shows results in the axisymmetric potential 
case. In this potential the Galactic inertial frame is employed 
to compute the orbits. For the 1642 stars, this figure shows 
the Lindblad diagram E versus /i, both constants of motion 
in the axisymmetric potential. Throughout this work, the 
units employed for h are 10^ km s~^ kpc, and 10^ km^ s“^ 
for E and J; the unit of distance is kpc. As it can be no¬ 
ticed in Figure [H in the axisymmetric case no structure is 
visible beyond an approximate homogeneous disc (toward 
the lower right part of the diagram), and an approximate 
homogeneous stellar halo (toward the centre-left part of the 
diagram). 

In a nonaxisymmetric potential, E and h computed 
in the Galactic inertial frame are not conserved along the 
stellar orbit; thus, these are not appropriate parameters 
to be compared in this case. After several combinations 
of orbital parameters, it was found that the particular 
combinations (F^min + Ema^)/^ versus (/imin + /lmax)/2 or 
(F^min + F^max)/2 vcrsus J gavc interesting diagrams. F^min, 
Fornax and /imin,/imax arc thc Corresponding minimum and 
maximum values of E and h along the stellar orbit. We 
call (F^min + Fornax)/2 thc ‘charactcristic’ orbital energy, and 
(hmin + /imax)/2 thc ‘characteristic’ orbital angular momen¬ 
tum. In Figure [2] the empirical diagram, characteristic en¬ 
ergy versus characteristic angular momentum, is presented, 
employing a nonaxisymmetric potential that includes the 
axisymmetric components, plus the central boxy bar (case 
(2) above). This plot shows a noticeable difference with re¬ 


spect to the axisymmetric case in Figure [T] In this case the 
observational points gather in several bands that represent 
the trapping regions produced by different resonances on the 
Galactic plane, as it will be shown in Section[5l Although our 
diagrams might seem similar to those of energy versus an- 
gular momentum employed bv lllelmi de Zeeuwl (|2000l i in 
their study of moving groups, the results and origin of the 
formed moving groups are fund amentally different. In the 
case of iHelmi de Zeeuwl (|2000l i. their method was applied 
to synthetic satellites created from the beginning into an ax¬ 
isymmetric simplified potential of the Galaxy. In this work, 
on the other hand, we produce a comprehensive study by ap¬ 
plying our method in a nonaxisymmetric Galactic potential, 
where some moving groups might be shaped by resonances 
generated by the Galactic bar. 

In Figure [3l an alternate form to the results in Figure 
is presented. This figure shows the characteristic orbital 
energy (F^min+F^max)/2 versus the orbital Jacobi constant J, 
which is conserved in the reference frame where the bar is at 
rest. In the same manner as in Figure [2l dots representing 
stars in the halo and disc agglomerate in noticeable lines 
or groups, conserving the same structure obtained in that 
figure. The crowded region of dots to the left in this diagram 
represents mainly disc stars, while the region to the right 
contains mainly halo stars. 

Figure S] shows how the diagram in Figure [2] changes 
when the spiral arms are included in the nonaxisymmetric 
potential. The difference between both diagrams is negli¬ 
gible for the case of halo stars, i.e., the spiral arms have 
little effect on halo stars, unlike the bar, where the effect 
on the formation of groups is clear, most likely due to its 
much larger mass (approximately ten times the mass of the 
spiral arms in the Milky Way case), and also because the 
perigalacticons of these stars approach close or within the 
bar region. However, the effect of the spiral arms is still of 
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Figure 1. The Lindblad diagram of orbital energy per unit 
mass versus the 2 ;-component of orbital angular momentum per 
unit mass for the 1642-star sample in the axisymmetric potential 
model of the Galaxy. The more bound (more negative orbital en¬ 
ergies) and largest positive z-component orbital angular momenta 
orbits, represented by the darkest region to the lower-right part of 
the plot, includes mainly Galactic disc stars, while the remainder 
are mostly halo stars. The units for h and E are 10^ km s“^ kpc 
and 10^ km^ s“^, respectively. 


Figure 2. Similar to Figure [T] here for the Galactic potential 
model including only the boxy bar in its nonaxisymmetric com¬ 
ponents. In this case the diagram plots a ‘characteristic’ orbital 
energy (in the Galactic inertial frame E is not conserved along a 
given orbit) versus a ‘characteristic’ orbital angular momentum 
{h is not conserved either), both given by calculating the average 
of the maximum and minimum values of E and h along a given 
orbit. In this case, noticeable lines of dots (agglomerations of stars 
in this diagram) are formed in the halo and disc regions, induced 
by the presence of the bar. 


importance on disc stars. In the following, results are given 
for the Galactic potential including only the bar as a nonax¬ 
isymmetric component; in a future study the effects of the 
spiral arms, in connection with the results presented in the 
next sections, will be analysed in detail. 

Figure [5] is a zoom of Figure [3l In this plot the agglom¬ 
erations formed in the disc (the left-ward crowded region) 
and in the halo (the remaining points) are illustrated in more 
detail. This figure is employed in the next section. 

The diagrams as in Figures [2] and [3] are the base of 
our method to delineate the agglomerations of points in the 
disc and halo regions. These agglomerations are associated 
with periodic orbits on the Galactic plane, as it is shown in 
the next Section [5l The next step is to investigate if known 
moving groups in the disc and halo are contained in these 
regions. In the following sections, some additional orbital 
properties are presented related with these diagrams. 


5 PERIODIC ORBITS ON THE GALACTIC 
PLANE 

To connect with the diagrams shown in the previous section, 
several families of periodic orbits on the Galactic plane have 
been computed in the nonaxisymmetric Galactic potential, 
that includes the axisymmetric components and the central 
boxy bar. These families of periodic orbits are generated by 
the presence of the bar, and are defined in the noninertial 
reference frame where the bar is at rest. In this frame, the 
x',y' plane is the Galactic plane, with the x' and y' axes 
pointing along the major and minor axes of the bar, respec- 



Figure 3. An alternate form to present the results in Figure 
[21 this diagram plots the characteristic orbital energy (Emin + 
Eniax)/2 versus the orbital Jacobi constant, in units of 10^ km^ 
s“^, computed in the reference frame of the bar. As in Figure (21 
noticeable agglomerations of stars, induced by the presence of the 
bar, are formed in the halo and disc regions. 
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(h, 


min~^^max 


)/2 



J 


Figure 4. As in Figure [2] but now including the spiral arms into 
the nonaxisymmetric barred potential. Both figures are similar in 
the case of the halo stars; i.e., the spiral arms have little effect 
on halo stars, contrary to the case for disc stars shown by the 
most crowded region to the right of the diagram, where some 
differences can be noticed. 



J 

Figure 5. This figure shows a zoom of Figure [3] to see in more 
detail the formation of the agglomerations in the disc and halo 
regions, and is at the same scale as the following Figure [G] for 
comparison. 

tively; the z' axis coincides with the z axis in the inertial 
reference frame. In this inertial frame the Galactic and bar 
rotations point in the clockwise direction. 

The use of Poincare-section diagrams has been useful 
to locate a starting periodic orbit in each family of peri- 
odic orbits, from which, using a Newton-Raphson method 
(|Press et ahlllQQ^ h we have computed how the orbit evolves 


Figure 6. Characteristic energy versus J for the nonaxisymmet¬ 
ric barred Milky Way model, showing in coloured solid lines and 
with Roman numbers from I to XI, different families of periodic 
orbits on the Galactic plane. These families are defined in the 
noninertial reference frame where the bar is at rest. The black 
points are the observational points plotted in Figure [5] Notice 
that several families coincide with crowded lines of dots in Figure 

m 

as the value of the orbital Jacobi constant is varied. For each 
computed orbit in a given family, the characteristic energy 
(^min + ^max)/2, cmploycd iu the previous section, and the 
orbital Jacobi constant J are obtained; thus, the entire fam¬ 
ily can be plotted in a diagram as in Figure O 

In Figure [6] several families of periodic orbits, marked 
with Roman numbers from I to XI, are plotted in a charac¬ 
teristic energy vs J diagram, along with the observational 
points in Figure [5] These families are shown by coloured 
solid lines. A given colour may represent more than one 
family. Notice that several of these families of periodic or¬ 
bits coincide with the crowded lines of dots (groups of stars) 
induced by the bar (see Figure [5|). This occurs even in the 
halo region of the background observational points, which is 
surprising, given that the periodic orbits lie on the Galactic 
plane. This influence of resonances extending with distance 
from the Galactic plane is analysed in Section [6l for specific 
families of periodic orbits in Figure [6l 

For stellar motion on the Galactic plane analysed in 
the noninertial reference frame, in Figure [7| we show the 
Poincare diagram in the plane x\vx' for the value of the 
Jacobi constant J= —2.3 units. Several resonant orbits are 
shown in this diagram. In Figures [8] and [9] we mark in this 
Poincare diagram the families of periodic orbits I to XI con¬ 
sidered in Figure [6l Notice that for this particular value of 
J the corresponding orbit in a given family can be stable or 
unstable, i.e. encircled or not by invariant curves. 

Figures [TO] to [20] show each individual family of peri¬ 
odic orbits given in Figure [6] In these figures the curve of 
the given family is plotted in two colours, red and green, 
which show the stable and unstable parts in each curve, re¬ 
spectively. The stability or instability of the computed peri- 
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Figure 7. For stellar motion on the Galactic plane, this figure Figure 9. Here we show the positions in the Poincare diagram 

shows the Poincare diagram in the plane a:', Ua,/for —2.3 units. for families VI to XI in Figure [6] Again, the colour of points 

corresponds to the colour of curves in Figure [6] 



X* (kpc) 


Figure 8. This is the left part of Poincare diagram in Figure [7| 
showing the positions of families I to V in Figure [6] The colour 
of points for each family is the same as that employed in Figure 
[6] for the corresponding curves. 

odic orb its in each fam ily was analysed following the method 
given byE enoni (|l965h . The agglomerations of points around 
a given family curve may occur in its stable parts. Also, in 
these figures five positions are taken in each family curve, 
labelled from 1 to 5, and their corresponding periodic orbits 
are shown in the five small frames, including the position 
number; the scale is the same in these frames, with distances 
in kpc. These orbits are plotted in the noninertial reference 
frame where the bar is at rest. 

As Figure [6] shows, all the family curves merge on both 


the left and right boundaries of the region where the points 
distribute. Figures [10] to [20] show that all the details of the 
periodic orbits tend to disappear at these boundaries, the 
orbits becoming almost circular. Due to the presence of the 
bar, this approximation to circular orbits becomes strictly 
true as the orbital characteristic energy increases. 

With respect to the inertial Galactic frame, the sense 
of orbital rotation (of the corresponding orbit in this frame, 
which in general will not be periodic) in these left and right 
boundaries in Figure [6] is as follows: in the right boundary, 
i.e., higher values of J, the orbital rotation is retrograde] 
in the left boundary, toward the region where the disc com¬ 
ponent distributes, this rotation is prograde. In Figures [TO] 
to [20] a blue starred point on each curve marks the approx¬ 
imate position where the sense of orbital rotation changes 
in the inertial frame, from retrograde to prograde and vice 
versa. This position corresponds approximately to the pe¬ 
riodic orbit in the family that crosses the Galactic centre 
(nearly localised by one of the five sample periodic orbits in 
the cases of families I,III,IV,V,VII,VIII,IX; see correspond¬ 
ing figures); this orbit will not have strictly h = 0, as is 
not conserved. 

An additional property of a given family of periodic or¬ 
bits, is that a point on its curve in a diagram as in Figure 
[6] i.e., characteristic energy versus J, is not necessarily as¬ 
sociated with a single periodic orbit. If a periodic orbit is 
not symmetric with respect to one of the axes of the bar, 
as in families III,VI,VIII,IX, the reflection of the orbit on 
this axis, which generates a different orbit in the given x',y' 
frame, will give the same point in the diagram as the origi¬ 
nal orbit; i.e., the characteristic energy and J are the same 
in both orbits. The existence of the reflected orbit is allowed 
due to the symmetries of the potential, i.e., the sense of both 
the x',y' axes can be changed simultaneously, recovering the 
original situation. Thus, both types of orbit coexist in the 
given x^,y^ frame. 
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Figure 10. This is family T’ in Fig. [6] The stable and unsta¬ 
ble parts along the curve are shown with red and green colours, 
respectively. Five periodic orbits along the family curve are illus¬ 
trated, numbered from 1 to 5 on the curve and in their frames. 
The orbits are plotted in the noninertial reference frame where 
the bar is at rest. The x' and y' axes on the Galactic plane point 
along the major and minor axes of the bar, respectively. The blue 
star on the family curve shows the approximate position where 
there is a change in orbital rotation as seen in the inertial Galac¬ 
tic frame: from retrograde to prograde going from right to left. 
In the inertial frame the Galactic and bar rotations point in the 
clockwise direction. Distances are in kpc. 






x' x’ x’ 

Figure 11. The same as Figure fTol here for the family labelled 
as TF in Figure [G] 


Figure 12. The same as Figure [TO] here for the family labelled 
as TIT in Figure [6] 



Figure 13. The same as Figure fTol here for the family labelled 
as TV’ in Figure [6] 


6 INFLUENCE OF RESONANCES ON 

Z-DISTANCE 

In Figure [6] some families of periodic orbits on the Galactic 
plane coincide with agglomerations of points in Figure (5] 
These agglomerations occur in the disc and halo regions; 
thus, the influence of these resonant families does not appear 
to be restricted necessarily to a region close to the Galactic 
plane. This influence can extend deep into the halo region. 

To show this result, in Figure!^ six agglomerations or 
regions, A,B,G,D,E,F in different colours, have been con¬ 
sidered from the observational points in Figure [5] Regions 
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Figure 14. The same as Figure fTOl here for the family labelled Figure 16. The same as Figure fTOl here for the family labelled 

as ‘V’ in Figure [G] as ‘VIF in Figure [G] 



x' x’ x’ 


Figure 15. The same as Figure fTol here for the family labelled 
as ‘Vr in Figure [G] 



x' x’ x’ 


Figure 17. The same as Figure [To] here for the family labelled 
as ‘VIIF in Figure [6] 


A,B,C are associated respectively with families IV,VII,VIII 
in Figure [6] (see Figures fT3lIT gIITT)) . Regions D,E,F are asso¬ 
ciated respectively with families V,IX,XI in Figure [6l which 
have stable parts in these regions (see Figures [TH [181 EQ!? 
but probably not with the corresponding neighbouring fam¬ 
ilies IV,VIII,X, which have unstable parts in these regions 
(see Figures [131 [13 fl^ • 

Figure [22l shows the distribution of the 1642-star sample 
on the plane of velocities U,V with respect to the inertial 
Galactic frame, U being negative toward the Galactic centre, 
and V positive in the direction of Galactic rotation. The 
points corresponding to stars in the regions A,B,G,D,E,F in 


Figure!^ are shown with their colours in that figure, along 
with a Bottlinger diagram showing some curves of constant 
orbital eccentricity, obtained with the approximation of an 
inver se square law for the radial force (|TrumDler V Weaverl 
Il953h . Regions D,E,F distribute in the prograde-retrograde 
transition middle zone, with high orbital eccentricities. 

For every point (star) in each of the regions 
A,B,G,D,E,F, the maximum orbital distance from the Galac¬ 
tic plane, |^|max, was computed, and corresponding diagrams 
1^1 max versus J were made for each region. These diagrams 
are shown in Figures [23lto[28l The z-distance can reach high 
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Figure 18. The same as Figure fTOl here for the family labelled Figure 20. The same as Figure fTOl here for the family labelled 

as ‘IX’ in Figure [G] as ‘XF in Figure [6] 



Figure 19. The same as Figure fTol here for the family labelled 
as ‘X’ in Figure [G] 

values, especially in the regions D,E,F for the halo compo¬ 
nent. 

For a sample of stars lying in regions D,E,F and reach¬ 
ing high ^-distances, Figures I29I30I31I respectively, show 
eight pairs of frames, four of them arranged in the first and 
second columns and the other four in the third and fourth 
columns. In each pair, the frame at the left gives the orbital 
projection on the Galactic plane x',y' in the reference frame 
where the bar is at rest, and the corresponding meridional 
orbit in the plane i?, z is shown in the frame at the right; 
R is the distance from the Galactic rotational axis z. As 
commented above, stars in these regions D,E,F are mainly 


Figure 21. Six regions A,B,C,D,E,F of observational points from 
Figure[5] shown in different colours. Regions A,B,C are associated 
respectively with families IV,VII,VIII in Figure[6](see Figures fT^ 
I1GII17|) . Regions D,E,F are associated with stable parts of families 
V,IX,XI in Figure [G] respectively (see F igures 11411181 1201) . 

associated with families V,IX,XI, respectively. This can be 
noticed comparing the orbital projections on the Galactic 
plane in these figures with the form of the periodic orbits in 
families V,IX,XI (see Figures fT^ fTSl [20l) . Figure [29] shows in 
the third pair of frames at the left side an orbit which seems 
to be associated with the family IV, and Figure [31] shows 
in the second and fourth pairs of frames at the left side 
two examples of orbits which seem to be associated with 
the family X, but their projections on the Galactic plane 
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Figure 22. The 1642-star sample in the U,V velocity plane. 
These velocities are with respect to the inertial Galactic frame; 
U < 0 toward the Galactic centre. Regions A,B,G,D,E,F in Fig¬ 
ure ED are shown with their corresponding colours. Some curves 
with constant orbital eccentricity in a Bottlinger diagram are 
also shown. From the inner to outer curves in the prograde and 
retrograde parts, the values of the eccentricity are 0.1,0.2,0.3, 
...,0.9,0.92,0.94,0.96,0.98 . 


have departed considerably from the unstable periodic orbit 
in this family X. The conclusion from Figures I29I30I311 is 
that some resonances on the Galactic plane can trap stars 
reaching high ^-distances from this plane. In t heir A /'-body 
models of barred galaxies, ICeverino fc Klvoinl (|2007h found 
that particles in the halo with a z-distance lower that 3 kpc 
could be trapped by the corotation and inner Lindblad res¬ 
onances. In our case, for the analysed resonances, the trap¬ 
ping z-distance can be larger. 

As it is shown in Figures [25] to 1281 for the regions 
C,D,E,F there is some internal separation of points for the 
1^1 max-distance. This separation is also present when we 
compute the orbits in the axisymmetric potential obtained 
by transforming the bar back into the original spherical 
bulge (see Section [3|. Figures [32] to [34] show the situation 
in the regions D,E,F. In these figures we plot |2:|max ver¬ 
sus Tmax, with Tmax the maximum orbital distance from the 
Galactic centre. The red points are obtained with the ax¬ 
isymmetric potential and the blue points with the nonax- 
isymmetric potential. 

The separation in |^|max in the regions C,D,E,F seems 
to be associated with an orbital transition between low- and 
high-|z| orbital domains. In Figure [35] we show how the max¬ 
imum |z|-amplitude varies in time for ten member stars in 
the region F; the present time is at t = 0. The red curves 
show the results with the axisymmetric potential and the 
blue curves the corresponding in the nonaxisymmetric po¬ 
tential. The first three rows in the figure show six orbits 
with points in the high-|z|max domain in the axisymmetric 
potential, but with points in the low-|2:|max domain in the 
nonaxisymmetric potential. Notice that for the first star in 


Figure 23. Maximum orbital distance from the Galactic plane 
for the stars in the region A, shown in Figure EH versus their 
corresponding J. 
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Figure 24. Same as in Figure [23] here for the region B of Figure 

El] 


the second row its orbit has suffered a transition to the low- 
1^1 max domain, but for the whole computed time of 10 Gyr it 
has been assigned to the high-jzjmax domain. The last two 
rows in the figure also show orbital transitions in the ax¬ 
isymmetric and nonaxisymmetric potentials. In particular 
for stars which can be trapped by the resonance associated 
to the region F, Figure [35] shows that the presence of the 
Galactic bar might favor the motion close to the Galactic 
plane, compared with that obtained in the axisymmetric 
potential; this result is also present in the resonance associ¬ 
ated to the region E, see Figure [33] The separations in |z|max 
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Figure 25. Same as in Figure[23l here for the region C of Figure Figure 27. Same as in Figure !^ here for the region E of Figure 


EH 


EH 


o 

a 



1^1111^11111p 


:—i -1-1-1-1-1-1-1-1-1-1-1- 


■ 

20 

_ ■ _ 

10 

Region D 


Region F 

8 

- 1 ■ - 
■ ■. ■ ■■ 

15 

_ ■ 


■ ■ 

" ■ 



6 

■ ■ 

L. - 

■ ■ ■ ■ ■ 

. ■ " > " ■ 

. ■ ■ 

o 

a 

^ 10 

■ 

— ■ — 



X 

CO 

■ ■ ■ H " 


■ ■ ■ - 

£ 


4 


N 

- 



5 

_ _ 

2 

— ■ ^ ^ — 


- 


> ■■■ ■■■■ >■ 


■ ■ _ 


■ ■■ ■ ■■ ' 



0 

_■ ■ _ 

_1__1_1_1__1_1_1__1_1_1_L 

0 

_ ■ ■ ■ _ 

1 1 1 1 1 1 1 1 1 1 1 1 1 


-1.8 -1.6 -1.4 -1.2 

-1.2 -1 -0.8 


J 


J 


Figure 26. Same as in Figure !^ here for the region D of Figure Figure 28. Same as in Figure!^ here for the region F of Figure 


EH 


EH 


might be explained by the existence of a major 1:1 resonant 
family of tube orbits which passes a few kpc above the Sun 
and which can be appreciated i n a number of the ‘verti¬ 
cal’ Poincare diagrams given in ISchuster Alleni (|l997l) . 
such as figs. 3a, 5a, and 6a. Such a family limits the phase 
space available for local halo-star orbits, and non-members 
of these tube orbits must navigate above or below the res¬ 
onant islands visible in these ‘vertical’ Poincare diagrams, 
leading to the separations visible in Figs. EH to EE and [32] 
to 1341 and in the behaviour of Fig. EH Also, stars which 
belong to such a 1:1 resonant family are not found in the so¬ 
lar neighbourhood, but further above or below the Galactic 


plane than the stars of the observational sample described in 
Section El maintaining these empty separations. This issue, 
and its possible consequences for the external shape of the 
Galactic halo, needs further study, analysing the resonant 
families external to the region F to see if they produce a 
similar behaviour. 

This section has shown that resonances on the Galactic 
plane can have a strong influence on stellar orbits reaching 
high-z distances from this plane. The stars of the Galac¬ 
tic halo may be strongly affected, acting as test particles 
influenced by massive nonaxisymmetric structures of the 
Galactic disc and bar. This effect may have important con- 
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Figure 29. For a sample of stars reaching high ^-distances and 
lying in region D in Figure EH this figure shows eight pairs of 
frames, four pairs are arranged in the first and second columns 
and the other four in the third and fourth columns. In each pair, 
the frame at the left gives the orbital projection on the Galactic 
plane x',y', in the reference frame where the bar is at rest, and 
the frame at the right gives the meridional orbit in the plane R, z. 



Figure 30. As in Figure [29l here for a sample of stars in region 
E in Figure [21] 


Figure 31. As in Figure [29] here for a sample of stars in region 
F in Figure [21] 
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Figure 32. For member stars in the region D, this figure shows 
in the axisymmetric potential (red points) and nonaxisymmetric 
potential (blue points) the maximum orbital 2 :-distance versus 
rmax, with rmax the maximum distance from the Galactic centre. 
The separation in |2;|max is present in both potentials. 



sequences in the formation of moving groups in the halo 
associated with periodic orbits on the Galactic plane. 

For stars with high-z distances, there are other fac¬ 
tors that may influence their stell ar kinematics, such as 
the triaxiality of the Galacti c halo (jPenarrubia et al]l2009l : 
iRoias-Nino et ahl 12012l . l2015l ) , that with the same resonant 
process as the nonaxisymmetric bar structures, trigger the 


production and support of stellar moving groups. The rela¬ 
tive importance of the kinematics imprinted by resonances 
of the nonaxisymmetric bar structure versus the one im¬ 
printed by the triaxiality of the halo is a study that will be 
presented in a future paper. 
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Figure 33. As in Figure [32l here for member stars in the region 

E. 
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Figure 34. As in Figure [3^ here for member stars in the region 

F. 
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7 SPECTRAL ANALYSIS 

To complement the analysis presented in previous sections, a 
spectral analysis of the three-dimensional orbits of the 1642- 
star sample has been made, and also for the two-dimensional 
orbits of the several families of periodic orbits on the Galac¬ 
tic plane compu ted in Section [5] We have considered the 
method given bv I Press et alJ (| 19921 ) for the spectral analysis 
of unevenly sampled data, and their fast subroutine fasper 
has been employed. In a given orbit, computed in the non- 
inertial reference frame where the bar is at rest, the spectral 
analysis was made for each Cartesian coordinate x'(t), y'(t), 
z'(t), or only x'(t), y'(t) in the case of the families of peri¬ 


Figure 35. Maximum orbital |z|-amplitude as a function of time 
for ten member stars in the region F. The red and blue curves 
show the results in the axisymmetric and nonaxisymmetric po¬ 
tentials, respectively. 


odic orbits. The dominant frequency (time“^) was obtained 
in each coordinate, i.e. the one with the largest amplitude 
or power in the spectrum. 

For the 1642-star sample. Figure [36] compares the ratios 
of the dominant angular fr equencies cc,/, uy (thi s figure 
is similar to figure 3.45 in iBinnev Tremaine! (|2008l )). The 
majority of the sample stars distributes around ujyi — 
1, or with the resonant condition cj-n = 0, with n=(l,—1,0). 
Other populated parts appear at the lower region in the fig¬ 
ure, especially at the lower left with n=(l,—3,0) and in the 
diagonal line n=(l,—1,—1). At the right side in this figure 
there is a concentration around n=(3,—1,0). The lower re¬ 
gion in this diagram corresponds to halo stars, in which cUz' 
decreases (the orbital period in the z-direction increases). 

To make a diagram similar to those constructed in Sec- 
tionUl in Figure [371 we plot for the 1642-star sample the ratio 
ujy! Iujx' versus the corresponding orbital Jacobi constant J 
(taking J ^ 1 as in Figure [5]). Again, several agglomerations 
of points are obtained, with particular details in the inter¬ 
val ~ (—2,—0.5) in J, where regions D,E,F in Figure!^ are 
located. 

The next step is to plot in a diagram like the one shown 
in Figure (371 the several families of periodic orbits computed 
in Section (5] Figures [38] and [39] show this diagram separately 
for each family; the black points in each frame correspond to 
the sample stars and the coloured lines to the given family. 
The number of the family is given in each frame and the 
employed colour is the same as in Figure [6] Notice that the 
ratio ujyf /ujx' changes in a given family as J is varied (when 
there is a transition between two levels in this ratio, its value 
may oscillate between the two levels and this appears as an 
overlapping of levels at the transition). The agglomerations 
of sample stars, which have three-dimensional orbits, are 
nearly covered by the families of periodic orbits, which lie 
on the Galactic plane. 
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Figure 36. Ratios of the dominant orbital angular frequencies 
in the three Cartesian coordinates x'(t),y'(t),z'(t) for the 1642 
sample stars. 


Figure 38. The ratio as a function of J for families 

of periodic orbits I to VI. The colour of lines in each family is 
that employed in Figure [6] The black points correspond to the 
1642-star sample. 



j 


Figure 37. The ratio uOy/ as a function of J for the 1642-star 
sample. Here we take J ^ 1 as in the zoom shown in Figure O 
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Figure 39. As in Figure !^ here for families VII to XL 


Now, in particular, the halo regions D,E,F shown in Fig¬ 
ure HD are considered. We have seen that these regions are 
likely associated with families V,IX,XI, respectively, which 
have stable parts in the positions of these regions. Thus, 
it is of interest to see what is the relation region-family: 
D-V, F-IX, and F-XI in a diagram Uy/jux' versus J. Fig¬ 
ures EOlllII and 1121 show the results. The black lines corre¬ 
spond to the family and the red points to the region. These 
figures show that there are approximate relations between 
the corresponding region and family. A three-dimensional 
orbit which may be trapped by a resonant family on the 


Galactic plane, will have a tube orbital projection on this 
plane which resembles the corresponding periodic orbit in 
the family. Only if this tube orbit is sufficiently narrow, the 
dominant frequencies in the x' and y' directions will be sim¬ 
ilar in both the projected and periodic orbits. The examples 
given in Figures [29l [30l and [31] show that some resulting 
tube orbital projections of three-dimensional orbits in re¬ 
gions D,E,F may have an appreciable thickness; thus, their 
dominat frequencies may differ from those obtained in the 
resonant family. This is shown explicitly in the next three fig¬ 
ures 43, 44, and 45 commented below. The conclusion from 
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J 


Figure 40. Comparison in a diagram uOy/ versus J between 
the family V of periodic orbits (black lines) and the region D in 
Figure [ 21 ] (red points). 



J 


Figure 42. Same as in Figure ^O] here for the family XI of peri¬ 
odic orbits and the region F in Figure [21] 



J 


SB and B21 is to plot directly the dominant frequencies f^^/, 
iy/ in the x' and y' coordinates versus J. This is shown in 
Figures SSI SSI and SSI The frequencies are given in (unit 
of time)“^ = (UT)“^, with 1 UT = 0.1 kpckm“^sec = 
9.784 X 10^ yr. In these figures the red and blue lines show 
respectively the values of ix' and f^/ for the given family of 
periodic orbits. The points (i.e. crosses) with magenta and 
cyan colours show the corresponding values of 4/ for stars 
in the given region. As commented above, and explicitly 
shown in these figures, ^x' foi* ^ three-dimensional orbit 
in a given region may not coincide with those obtained in 
the associated resonant family. A more complete comparison 
would be obtained with these type of figures by plotting not 
only the dominant frequencies in each case, but also other 
peak frequencies in the corresponding spectrum, especially 
if there are peaks with power similar to that obtained in the 
dominant frequency and thus not considered in the previous 
figures. 


8 APPLICATION OF THE METHOD 


Figure 41. Same as in Figure l40] here for the family IX of peri¬ 
odic orbits and the region E in Figure [21] 


this analysis of frequencies is that a diagram ujy/ jux' versus 
J may not give a complete answer concerning the trapping 
of a three-dimensional orbit by a resonant family on the 
Galactic plane; it is necessary to complement the analysis 
with figures showing the orbital projection on this plane, and 
include other frequencies in the spectrum, as commented in 
the following. As concluded in the previous section, and com¬ 
plemented with the analysis in this section, the halo regions 
D,E,F have stars that are trapped by resonant families on 
the Galactic plane. 

An alternate form to present the results in Figures [40] 


To illustrate the possible relation between the agglomera¬ 
tions of stars trapped by resonances obtained in the empir¬ 
ical diagram of characteristic energy versus J presented in 
Section 2] and groups of stars catalogued as moving groups, 
the position of the Kapteyn group, a known moving group 
in the Galactic halo, is considered for this diagram. The 
1642 halo stars and disc stars in the solar neighbourhood 
presented in Section [2] are again employed to delineate this 
diagram; they serve as a reference to compare with the cor¬ 
responding positions in the diagram of the sample stars in 
the Kapteyn group. In fact, some of these 1642 stars are 
members of this moving group. 

Eighteen member stars of the Kapteyn moving gro up 
hav e been selected from Eggen I (|l965al . Il977l . Il990l . ll996al lH) 
and IWvlie-de Boer et al.l (|2Q10l ). The positions, proper mo- 


© 0000 RAS, MNRAS 000 , 000-000 




















Resonant Trapping in the Galactie Dise and Halo and its Relation with Moving Groups 17 



J 

Figure 43. The frequencies in (UT)“^, with 1 UT = 0.1 

kpckm“^ sec, for the family V of periodic orbits and the region D 
in Figure [21] The red and blue lines show respectively the values 
of ix' i^y' for fho family V and the points with magenta and cyan 
colours show respectively the values of fa;/,f^/ for the region D. 



Figure 44. As in Figure 03] here for the family IX and the region 

E. 

tions, radial velocities, and distances of these stars have been 
improved as was performed for the sample of 1642 stars. 

Figure liHl shows the diagram with the observational 
points employed in Section 0] and some curves of families 
of periodic orbits of Section (5] along with red points corre¬ 
sponding to stellar members of the Kapteyn moving group. 
Comparing with Figure (6] these members arrange mainly 
around families III, IV, and V. However, only the family V 
has a stable part in its corresponding region covered by these 



J 

Figure 45. As in Figure03] here for the family XI and the region 

F. 



Figure 46. Diagram of characteristic energy versus J in the non- 
axisymmetric Galactic potential with the bar, showing with red 
points the positions of 18 members of the Kapteyn moving group. 
The small black dots correspond to the sample stars employed in 
Section |4] the black curves are some families of periodic orbits 
given in Section [5] The orbits of points (stars) marked as 1 to 6 
are presented in Figures [47] and 1481 

member stars; thus, the points marked as 1 to 6 in Figure [46] 
might have at the present time a relation with this resonant 
family V. 

To see if this relation exists, for the six points (stars) 
marked as 1 to 6 in Figure [46] members of the Kapteyn mov¬ 
ing group. Figures 07] and 08] show at the left side their corre¬ 
sponding orbital projections on the Galactic plane, plotted 
in the x^,y' frame where the bar is at rest. Their meridional 
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Figure 47. For the three points marked as 1,2,3 in Figure H6l the 
frames on the left show the corresponding orbital projections on 
the Galactic plane in the x',y' axes, and the frames on the right 
show the meridional orbits. 


orbits are shown on the right side. Points 1,2,4 have an or¬ 
bital projection on the Galactic plane which is a tube orbit, 
of the type that would be obtained around stable periodic 
orbits in the family V (see panel 2 of Figure [M]), in accor¬ 
dance with their positions on this family, as shown in Figure 
M On the other hand, points 3 and 6 are very close to the 
curve of the family V, but their orbital projections on the 
Galactic plane resemble an orbit departing from the unsta¬ 
ble periodic orbits in the family IV (see frame 3 in Figure 
[IS]). The orbit corresponding to point 5 does not seem to 
be associated with any resonance. Thus, the majority of the 
member stars of the Kapteyn moving group does not appear 
to have a relation with a resonant family on the Galactic 
plane. 

On the other hand, a direct connection between the 
Kapteyn moving group and the well studied globular cluster 
u Gentauri has been argued bv lWvlie-de Boer et al.l (|2010l l 
using detailed kinematical and chemical analyses. They pro¬ 
pose that many, or all, members of this moving group have 
been left behind as remnants as this globular cluster (ex¬ 
dwarf g alaxy) was dragge d into the Galaxy by dynamical 
friction. iMeza et al.l (|2005l ) have carried out numerical sim¬ 
ulations of such a process, and find a local group of stars 
which seem to relate to uj Gen both in their angul ar momenta 
and in their chemical abundance characteristics. ISilva et"aI1 
(|2012l ) use the Bottlinger, Toomre and [Fe/H],V(rot) dia¬ 
grams to identify in the solar neighbourhood halo moving- 
group candid ates relating to both t he Kapteyn and the cu 
Gen entities. ISchuster et al. I (|2012l ) conclude that u Gen 
may have contributed in a significant way to their ‘low- 
alpha’ local halo component. If the Kapteyn group is of 
‘cosmological’ origin, i.e. all stars having formed together 
in uj Gen, then the appearance of Figure 46 may be due to 
the migration of the Kapteyn stars to the nearby resonant 
families III, IV, and V during their trajectories while being 


Figure 48. As in Figure ITfl here for the three points marked as 
4,5,6 in Figure |46] 

lost to the Galaxy from lu Gen. That is. Figure 46 seems 
to suggest a possible secular origin for the Kapteyn mov¬ 
ing group, i.e. a product of the resonant trapping by the 
Galactic bar, but this is misleading, as shown here. 

In a second paper, other moving groups in the Galactic 
halo will be analysed using the diagram of characteristic 
energy versus J. The method will also be applied to some 
known moving groups in the the Galactic disc. 


9 CONCLUSIONS 

In this study we have analysed the trapping of stars by res¬ 
onances on the Galactic plane created by the Galactic bar. 
Our aim is to relate this mechanism with moving groups in 
our Galaxy, especially with moving groups in the Galactic 
halo. A new method is presented to delineate the trapping 
regions. To show how the method works, a sample of halo 
stars and disc stars in the solar neighbourhood has been 
employed. The orbits of these sample stars are computed 
in a nonaxisymmetric Galactic potential, and empirical di¬ 
agrams are constructed using some particular orbital prop¬ 
erties. These diagrams plot a characteristic energy versus a 
characteristic angular momentum, or a characteristic energy 
versus the orbital Jacobi constant, if the Galactic potential 
employs only the bar as its nonaxisymmetric component. 

Some agglomerations of points in the disc and halo re¬ 
gions are obtained in these diagrams, which coincide with 
some families of periodic orbits on the Galactic plane, in¬ 
duced by the Galactic bar and plotted in the same diagrams. 
The influence of these resonant families can extend some kpc 
from the Galactic plane, explaining the observed coincidence 
in the halo region. Our future aim with this analysis is to 
investigate if these agglomerations obtained in the disc and 
halo can be the sites where some known moving groups lie. 

The proposed method has been illustrated here with the 
Kapteyn group, a known moving group in the Galactic halo. 
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resulting in that a majority of member stars in this group 
does not have a relation with any of the resonant families 
or periodic orbits obtained in the employed Galactic poten¬ 
tial, instead, we propose that these group might rather have 
an origin on the dissolution of a dwarf galaxy or a globu¬ 
lar cluster. In a second paper the method will be applied 
to other groups in the Galactic halo, including also some 
known moving groups in the disc, to provide one piece of 
information on their possible, or not, resonant origin. 
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